Statistical mechanics of bend flexoelectricity and the twist-bend phase in bent-core 

hquid crystals 
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We develop a Landau theory for bend flexoelectricity in liquid crystals of bent-core molecules. 
In the nematic phase of the model, the bend fiexoelectric coefficient increases as we reduce the 
temperature toward the nematic to polar phase transition. At this critical point, there is a second 
order transition from high-temperature uniform nematic phase to low-temperature nonuniform polar 
phase composed of twist-bend or splay-bend deformations. To test the predictions of Landau theory, 
we perform Monte Carlo simulations to find the director and polarization configurations as functions 
of temperature, applied electric field, and interaction parameters. 



I. INTRODUCTION 

Bent-core liquid crystals exhibit a rich variety of phases 
with different types of orientational order. At high tem- 
perature they form an isotropic phase, with disorder in 
all of the molecular axes. As the temperature decreases 
they form a nematic phase, in which the long axes of 
the molecules tend to align up or down along a director 
fi. The transverse orientations of the molecules can then 
have various types of order in the plane perpendicular to 
n: they might be disordered (leading to a uniaxial ne- 
matic phase), or have nematic order (leading to a biaxial 
nematic phase) , or have polar order (leadirig to a net po- 
larization). A classic argument of Meyer [l| shows that 
polar order of the transverse directions couples to bend 
variations in the main director fi. As a result, it is par- 
ticularly easy to induce polar order and director bend in 
bent-core liquid crystals, compared with analogous rod- 
like liquid crystals. 

Recent research has found two remarkable physical 
phenomena arising from polar order and director bend 
in bent-core liquid crystals: 

a. Flexoelectricity: The fiexoelectric effect is a linear 
coupling between polar order and director deformations 
in the uniaxial or biaxial nematic phase. An imposed 
deformation of the nematic director leads to an electro- 
static polarization, and conversely, an applied electric 
field induces a deformation of the nematic director. Since 
the initial discovery by Meyer this effect has drawn 
great interest due to its possible application in devices 
such as display panels, actuators, micropower genera- 
tors, and electro-mechanical transducers for sensing or 
energy-harvesting uses. In general, the director deforma- 
tion might be either splay or bend, but one would expect 
the bend fiexoelectric effect to be dominant in bent-core 
liquid crystals. Indeed, recent experiments by Harden et 
al. 0, 0] have found that bent-core liquid crystals have a 
surprisingly large bend fiexoelectric coefficient, up to 35 
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nC/m, roughly three orders of magnitude larger than the 
typical value of 3-20 pC/m in rod- like liquid crystals. 

For flexoelectricity, the key theoretical question is how 
to explain the large effect found in bent-core liquid crys- 
tals, so that it can be exploited for technological appli- 
cation. In a previous paper we conjectured that the 
large fiexoelectric effect is a statistical phenomenon as- 
sociated with nearby polar phase. Near a polar phase, a 
nematic liquid crystal is on the verge of developing spon- 
taneous polar order, and hence any deformation of the 
director should induce a large polar response. In the pre- 
vious paper, we explored that concept by constructing a 
model for splay fiexoelectricity in a system of pear-shaped 
molecules. It is still necessary to extend the model to the 
more complex but physically realistic case of bend flexo- 
electricity in bent-core liquid crystals. 

We recognize that experimental measurements of fiex- 
oelectricity in bent-core liquid crystals are controver- 
sial d, [1]. Regardless of the experimental fiexoelectric 
coefficient in any particular material, we would like to 
understand what behavior is theoretically possible. 

b. Modulated phases: In an influential theoretical 
paper Q , Dozov noted that bent-core liquid crystals have 
a very low energy cost for bend in the director, and de- 
scribed this low energy cost by a very small bend elastic 
constant K^. He then speculated that K3 could actu- 
ally become negative in some bent-core liquid crystals. 
(In this case, the free energy would need to be stabilized 
by higher-order terms in the bend.) Dozov showed the- 
oretically that the uniform nematic phase would become 
unstable to the formation of a modulated phase, which 
could be either a twist-bend or a splay-bend phase. Re- 
lated simulations have been done by Memmer [8]. This 
theoretical prediction is now supported experimentally 
by optical observations of spontaireous periodic deforma- 
tions in bimesogens 0, and by extensive studies of a 
liquid crystal dimer, which identify the twist-bend phase 
using a series of techniques including small-angle x-ray 
scattering, modulated differential scanning calorimetry, 
dielectric spectroscopy, and magnetic resonance [lOj . 

For the twist-bend and splay-bend modulated phases, 
there are two important theoretical questions. First, how 
does the elastic constant become negative; what is 
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the physical meaning of this negative value? Second, 
what happens to the polar order in these phases? The 
schematic illustrations in Dozov's paper Q clearly show 
the presence of local polar order, but this polar order is 
not included in his theoretical formalism. 

The purpose of this paper is to develop a theory of 
polar order and director bend in bent-core liquid crys- 
tals, which can explain both flexoelectricity and modu- 
lated phases in these materials. For this work, we use 
both Landau theory and lattice simulations. In Landau 
theory, we construct a free energy functional in terms 
of director gradients and polarization. In the uniform 
nematic phase, we can minimize the free energy to find 
the optimal polarization for fixed bend, or conversely the 
optimal bend for fixed polarization (or applied electric 
field). This minimization gives the bend flexoelectric co- 
efficient 63, as well as the effective bend elastic constant 
and effective dielectric coefficient. In particular, it shows 
the difference between the bare elastic constant and 
the effective (or renormalized) elastic constant Kf^. As 
the temperature decreases toward a critical temperature, 
63 increases and K^^ goes to zero. Below that critical 
temperature, the uniform nematic phase becomes unsta- 
ble to the formation of a new phase which has both po- 
lar order and spontaneous bend, and hence must have 
a modulated structure. By minimizing the free energy 
over variational forms for the director and the polariza- 
tion, we find that the director configuration is equivalent 
to Dozov's twist-bend or splay-bend phase, and we also 
determine the accompanying polarization configuration. 

In simulations, we construct a lattice Hamiltqnian that 
generalizes the classic Lebwohl-Lasher model [H,] for ne- 
matic liquid crystals by including two vectors on each 
site, which represent the long molecular axis and the 
transverse direction, respectively. This model is analo- 
gous to our previous lattice model for splay flexoelectric- 
ity [3 , but now extended to bend flexoelectricity. We run 
Monte Carlo simulations of this Hamiltonian in both the 
high-temperature uniform nematic phase and the low- 
temperature modulated polar phase. These simulations 
give results that are equivalent to Landau theory, but 
without the mean-field and variational assumptions of 
Landau theory. In the uniform nematic phase, they show 
the increasing flexoelectric coefficient as the temperature 
decreases. In the low temperature phase, they show the 
structure of the modulated phase, which may be twist- 
bend or splay-bend, depending on model parameters. 

The plan of this paper is as follows. In Sees. II and III, 
we develop a theory for the ffexoelectric effect in the uni- 
form nematic phase, first using Landau theory and then 
lattice simulations. In Sees. IV and V, we extend the the- 
ory to the low-temperature modulated polar phase, now 
in the reverse sequence of lattice simulations and then 
Landau theory. In Sec. VI we discuss and summarize 
the conclusions of this study. Finally, in the Appendix 
we show how the lattice model can also be solved using 
mean-field theory. 



II. FLEXOELECTRIC EFFECT: LANDAU 
THEORY 

To describe the flexoelectric effect in a uniform nematic 
phase, we must construct the free energy density in terms 
of the nematic director n and the polarization P. We 
assume that gradients of n are small, and P is also small, 
so that we can work to quadratic order in both of these 
quantities. (We will go to higher order in a later section, 
when discussing modulated phases.) The free energy has 
three parts: 

a. The free energy of director gradients is just the stan- 
dard Oseen-Frank free energy. It can be written as 

F„„ = ixiS2 + l/^aT^ + \k^^\ (1) 

in terms of the splay vector S ~ n(V ■ n), the twist 
pseudoscalar 7" = n • (V x n), and the bend vector 
B = n X (V X n). 

b. Polar order does not occur in the uniform nematic 
phase, and hence it must have a positive free energy cost. 
To lowest order, this cost can be written as 

Fpp = i^^P", (2) 

where /i is an arbitrary coefficient. There may be sev- 
eral physical contributions to /i. Entropy makes a posi- 
tive contribution proportional to temperature T, and the 
electrostatic energy makes another positive contribution 
independent of T . Because we are modeling bent-core 
liquid crystals, we suppose that there are packing con- 
siderations that favor polar order, in competition with 
entropy and electrostatics. This packing energy can be 
modeled by a negative part of the free energy, propor- 
tional to P^, independent of T. Hence, the combination 
of these effects gives a coefficient /x that varies with tem- 
perature, and can be written as /x = /Li'(T — Tq). Note 
that T — Tq is positive in the uniform nematic phase. 

c. Following the argument of Meyer polar order is 
coupled to both splay and bend in the director. Wc al- 
ready considered the splay coupling in our previous paper 
on splay fiexoelectricity \4^. In this paper we are con- 
cerned with bend flexoelectricity, and hence we consider 
the coupling 

F„p = -AB • P, (3) 

which favors polar order along the bend direction, per- 
pendicular to the director, and hence describes ordering 
of the transverse orientations of the molecules. 

Putting these three pieces together, we obtain the total 
free energy density 

F = ^KiS^ + ^K^T^ + iifgB^ + ^A^P^ - AB • P (4) 

to quadratic order in director gradients and polar order. 
Note that the last three terms are a quadratic form in 
B and P. The uniform nematic phase is only stable if 
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the quadratic form is positive-definite, whicfi occurs if 
< /ii^a- This condition can be rewritten as 



^3 



or equivalently as 



T>Tr.=Tn 



(5) 



(6) 



Below that critical temperature, the uniform nematic 
phase must become unstable to a phase with director gra- 
dients and polar order, which will be discussed in Sees. IV 
and V below. 

To model the flexoelectric effect, we minimize the free 
energy of Eq. Q over the polarization P for fixed bend 
B to obtain 



P = egB, 

where the bend flexoelectric coefficient is given by 

A 



A 

63 = - = 



(7) 



(8) 



This coefficient increases toward a finite limit as the tem- 
perature decreases toward the critical temperature Tc- 
Substituting the expression ([7]) for P back into Eq. 
gives the effective free energy 
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(9) 



where the effective, renormalized bend elastic constant is 
given by 



K. 



eff 



A^ 



A2 



A*' (7^ - To) 



(10) 



Note that there is an important physical distinction be- 
tween the bare elastic coefficient K3 and the renormalized 
coefficient Kf^: the bare coefficient K3 gives the energy 
cost of a bend in a hypothetical experiment where the 
polarization is constrained to be zero, while the renor- 
malized coefficient Kf^ gives the energy cost of a bend 
in an experiment where the polarization is free to relax to 
its lowest-free-energy value. This distinction is essentially 
the same as the distinction between the bend coefficient 
at constant electric displacement D compared with 
at constant electric field E, as discussed by Castles 
et al. in the context of blue phase stability ^12]. The 
renormalized coefficient is apparently the coefficient K^^ 
calculated by Cestari et al. through molecular field the- 
ory with atomistic modeling, because their calculation 
involves averaging over molecular distributions that re- 
spond to bend [ij]. The bare coefficient is always 
positive, while the renormalized coefficient K"^^ passes 
through zero as the temperature passes through Tc- In 
Sec. V below, we will discuss the behavior below Tr. 



and will show that K^^ is the negative bend constant 
discussed by Dozov [7|] . 

To model the converse flexoelectric effect, we add an 
electric field coupling to the polarization into the free 
energy ([4]), which gives 

= ixiS2 + lifaT" + i/^3B2 + i/ip2 - AB • P - E • P. 

(11) 

We minimize this free energy over both bend B and po- 
larization P for fixed electric field E to obtain 



B 
P 



AE 



Mi^a - A2 ■ 
fiKs - A2 ■ 



(12a) 
(12b) 



Equation (jl2al) shows that an electric field induces a bend 
through the converse flexoelectric effect, while Eq. (|12bp 
shows that an electric field induces a polarization. The 
latter equation can be compared with the standard ex- 
pression for the induced polarization in a dielectric ma- 
terial, F = eo{e — 1)E. This comparison shows that the 
effective dielectric constant is 



,cff 



1 



1 



eo(A*-AV^3) 



(13) 



Hence, the effective dielectric constant is renormalized 
upward by the coupling between polarization and bend, 
compared with the value e = 1-f l/(eoM) without this cou- 
pling. At the critical temperature Tc where the uniform 
nematic phase becomes unstable, the effective dielectric 
constant diverges. 

Here, we should emphasize which quantities diverge 
and which remain finite at the critical point T^. When 
we apply an electric field in the converse fiexoelectric ef- 
fect, the induced bend and polarization are given by sus- 
ceptibilities multiplying the field. These susceptibilities 
both diverge as T —> Tc, as usual for susceptibilities at 
a second-order phase transition. Likewise, if we were to 
apply a torque that couples directly to the bend, the in- 
duced bend and polarization would be given by suscepti- 
bilities multiplying the torque, and those susceptibilities 
would also diverge as T — However, the flexoelec- 
tric coefficient 63 is not exactly a susceptibility; rather, it 
is the ratio between bend and polarization, which both 
diverge at the transition. This ratio does not diverge, 
but approaches the finite maximum value ef^^ — K3/X. 
By comparing this maximum value with the bend elastic 
constant and dielectric constant as T ^ Tc, we see that 



= X3eo(e-l) = /f3°%(e'=«-l). 



(14) 



This equation for the maximum flexoelectric coefficient 
is approximately the same as the limit derived by Castles 
et al. However, they derive the limit based on argu- 
ments about the conservation of energy. We would not 
say that it is related to energy conservation, but rather 
that it is the limit of stability of the uniform nematic 
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FIG. 1: Lattice model for neighboring bent-core liquid crystal 
molecules, showing possible orientations of two neighboring 
molecules under a bend deformation. 



phase. Beyond that point, the uniform nematic phase 
becomes unstable to the formation of a modulated phase, 
as will be discussed in later sections. 



III. FLEXOELECTRIC EFFECT: LATTICE 
MODEL 

In the previous section, we developed a Landau theory 
for the bend flexoelectric effect near a transition to an in- 
cipient polar phase. In this section, we develop a lattice 
simulation model to describe the same effect. This model 
is an extension of the lattice model for splay flexoelec- 
tricity in our previous paper [4] . The lattice simulations 
will provide clear visualizations of the types of molecu- 
lar order that occur in the flexoelectric effect. They will 
also allow us to avoid the standard limitations of Landau 
theory, which neglects correlated fluctuations and applies 
only over a limited temperature interval. 

We consider a simple cubic lattice with a bent-core 
liquid crystal molecule on each site. The orientation of 
each molecule is characterized by two orthogonal unit 
vectors, as shown in Fig. [TJ The vector hi represents 
the long molecular axis, while the vector hi represents 
the transverse dipole of the molecule at site i. In the 
uniaxial nematic phase without any bend, the n vectors 
are ordered but the b vectors are disordered. We would 
like to model how a bend in n induces polar order in 6, 
or conversely, how an electric field applied to b induces 
bend in n. 

For the lattice interaction between molecules on neigh- 
boring sites i and j, we must consider several terms. 
First, we need the standard Lebwohl-Lasher [11] inter- 
action —A{ni ■ fij)'^, which favors nematic order of the 
molecules. Second, to describe a bent-core liquid crys- 
tal with incipient polar order, we need an interaction of 
the form —Bibi ■ bj. This term favors polar order of the 
transverse dipoles, driven by packing energy or any other 
mechanism. Third, there may also be a higher-order in- 
teraction of the form —B2{bi ■ bj)'^, which favors nematic 
order of the transverse dipoles, perpendicular to the main 



nematic axis. This term gives the possibilit y o f a biaxial 
nematic phase, as in the model of Straley Fourth, 
there must be an interaction ■ bi between each dipole 
and the applied electric field. 

Finally, we need a coupling of b with the local bend in 
h. Following the argument of Ref. [1], a lattice version 
of the bend vector between sites i and j can be written 
as 



[fi X (V X n) 



1 r 
2 



+hi{hi ■ nj){fij ■ hj) (15) 
-hj{rii ■ fij){rij ■ rii 



where — {rj —i'i)/\rj — r.i \ is the unit vector from site 
i to j on the lattice. Note that this expression is invariant 
under the symmetry operations i <r^ j, fii ^ —fii, and 
hj — > —fij. This expression for the local bend can be 
coupled with the local polar order, averaged over sites i 
and j, to give 

Vbcnd = C[n X (V X n)]jj ■ ^ijtA 
C r - 

-{b^ ■ rij)[r^j ■ {hj + n^n, • hj)}^ . (16) 

Combining all the terms, our final expression for the 
lattice Hamiltonian is 



H = - 



E 



A{h, ■ hjf + Bib, ■ bj + B2{b, • bjf 
C r - 

- |^(6j • hi)[r,j ■ {h, + hj{h, ■ h^)}] 
-{bi ■ hj)[fij ■ {hj + h,{hi ■ hj)}] 



(17) 



This expression is analogous to the lattice Hamiltonian 
in Ref. [1] , except that it has two orthogonal unit vectors 
on each site instead of just one, so that it can describe 
bend instead of splay flexoelectricity. 

We carry out Monte Carlo simulations of a system of 
bent-core molecules interacting with the Hamiltonian of 
Eq. (|17p on a simple cubic lattice of size 16 x 16 x 16. 
When an electric field is applied, it is in the Z-direction. 
The simulation box has periodic boundary conditions in 
Z , but the boundaries in X and Y are free so that the 
system can form bend across those directions. In each 
Monte Carlo step, a molecule is chosen randomly and it 
is slightly rotated about a random axis. The change in 
energy /S.E is calculated, and the step is accepted or re- 
jected following the Metropolis algorithm. Starting from 
a high-temperature isotropic state, the system is cooled 
down slowly with temperature steps of AT = 0.01. The 
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(a) 





FIG. 2: Equilibrium configurations from Monte Carlo simula- 
tions showing (a) isotropic phase, (b) uniaxial nematic phase, 
(c) biaxial nematic phase, and (d) polar phase. Cross sec- 
tions through the XZ-plsaie are shown, and the color on each 
molecule represents the orientation of its short axis b with 
respect to the laboratory Z-axis. 



final configuration at each temperature is taken as the ini- 
tial configuration for the next lower temperature. Typ- 
ical runs take about 10^ steps to come to equilibrium, 
while runs near phase transitions take about 6 x 10^ steps. 

Figure [2] shows snapshots of the molecular configura- 
tion in four phases at zero electric field. At high tem- 
perature, the system is in an isotropic (/) phase, with 
disorder in both n and b. As the temperature decreases, 
it forms a uniaxial nematic {Nu) phase, with nematic 
order in n but disorder in b vectors, which are uniformly 
distributed in the plane perpendicular to the average di- 
rector. At lower temperature, it forms a biaxial nematic 
{Nb) phase, with nematic order in both h and b. In the 
Nb phase, the b vectors have two favored orientations 
perpendicular to the average director. At the lowest tem- 
perature, it forms a polar (P) phase, with nematic order 
in n and polar order in b vectors, which now have one fa- 
vored orientation perpendicular to the average director. 
In this case, the the polar order induces bend in the di- 
rector. (At longer length scales more complex modulated 
structures are seen, as discussed in the following section.) 

To characterize the orientational order of these phases, 
we use a standard set of four uniaxial and biaxial nematic 
order parameters, supplemented by an additional polar 
order parameter. The nematic order is represented by 
the supermatrix [l5l - [l7| 



qAB _ 
'-'ab — 




0.0 0.2 0.4 0.6 0.8 1.0 

temperature(T) 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

temperature(T) 



FIG. 3: Monte Carlo simulation results for the order pa- 
rameters S (>), V (x) and P (o) as functions of tempera- 
ture T, for difi'erent values of the coupling coefficients cho- 
sen to show different types of transitions, (a) Bi — 0.04, 
B2 = 0.3. (b) Bi = 0.12, B2 = 0.3. (c) Bi = 0.4, B2 = 0.3. 
(d) Bi = 0.04, B2 = 0.95. In all cases A = 1.0, C = 0.3, and 
E = 0. 



Here, the subscripts a and b denote molecular axes while 
A and B represent laboratory axes, and laA is the di- 
rection cosine between molecular and laboratory axes, 
laA = a ■ A. In the principal coordinate system of this 
supermatrix, the nematic order parameters can be ex- 
pressed as |15l-[l7| 



S 
V 



qZZ 

'-'zz 

1 

3 



:;ZZ 



T=S 



XX 

zz 



S 



YY 



:;_\(cXX _ cYY\ _ (qXX _ cYY\i 
o [\ XX ^xx ) y yy yy -'J 



(19) 



(18) 



As a physical interpretation, S measures the ordering of 
the molecular z-axis with respect to the laboratory Z- 
axis and is the usual nematic order parameter, while U 
measures the difference in ordering of the molecular x and 
y-axes with respect to the laboratory Z-axis. Conversely, 
T measures the difference in ordering of the molecular 
z-axis with respect to the laboratory X and F-axes. Fi- 
nally, V measures the ordering of the molecular x and 
y-axes with respect to the laboratory X and F-axes. Of 
these four parameters, S and U are non-zero in the uni- 
axial and biaxial phases, while T and V vanish in the 
uniaxial phase but are non-zero in the biaxial phase. In 
the limit of low temperature and high orientational or- 
der, S and V tend to 1, while U and T tend to 0. We 
calculate all four of these order parameters following the 
procedure of Refs. (T5l - [T7| , but we only present results 
for the uniaxial nematic order parameter S and biaxial 
nematic order parameter V. 

In addition to the nematic order parameters, we define 
the polar order parameter as the average of the molecular 
dipole vectors, P = (b), and we present results for the 
magnitude P = \P\. 

The order parameters 5, and P are time-averaged 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 



temperature(T) 



temperature(T) 
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(d) 













0.0 0.2 0.4 0.6 O.S 1.0 1.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 
temperature(T) temperature(T) 



FIG. 4: Monte Carlo simulation results for the bend per bond 
as a function of temperature T, for different values of the 
coupling coefficients, (a) Bi = 0.04, B2 = 0.3. (b) Bi = 0.12, 
B2 = 0.3. (c) Bi = 0.4, B2 = 0.3. (d) Bi = 0.04, B2 = 0.95. 
In aU cases ^ = 1.0, C = 0.3, and E ^ 0. 



in the production cycle of simulations. Plots of the or- 
der parameters for different values of the coupling coeffi- 
cients are shown in Fig. [31 In all these cases, the isotropic 
phase is stable at high temperatures, followed by phases 
of reduced symmetry at low temperatures. We find four 
types of transitions depending upon the relative coupling 
coefficients: (i) For small biaxial coupling B2 and even 
smaller polar coupling Bi , there is a series of transitions 
/ — > Njj —5- Nb —5- P- At each transition, the system 
acquires additional order, (ii) For somewhat larger polar 
coupling Bi, a direct transition from uniaxial nematic 
to polar phase occurs, with no intermediate biaxial ne- 
matic phase, / — i> Njj P. (iii) For even stronger Bi, 
a direct transition from isotropic to polar phase takes 
place, with no intermediate uniaxial or biaxial nematic 
phase, I ^ P. (iv) If the molecules have a very strong 
biaxial coupling B2, the transition occurs goes directly 
from isotropic to biaxial nematic, with no intermediate 
uniaxial nematic phase, / — >■ Nb P- From these ob- 
servations, we see that the stability of the biaxial nematic 
phase can be tuned with the polar strength Bi and bi- 
axiality B2 of the molecules. 

A key feature of our model is that bend in n is coupled 
to polar order in b. For that reason, we expect nonzero 
bend to occur spontaneously as the system moves into 
the polar phase, even under zero applied electric field. 
To test this concept, we calculate the average bend from 
Eq. (1151) as a function of temperature and coupling coef- 
ficients. The results are shown in Fig. S) Note that this 
spontaneous bend is zero in the isotropic, uniaxial ne- 
matic, and biaxial nematic phases, and becomes nonzero 
in the polar phase. One further point about the order 
parameters should be noted here: We calculate S, V, 






Isotropic 










(b) 


Polar 



0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 5: Phase diagram of the model system as a function of 
(a) T and B2, for fixed Bi = 0.03. (b) T and Bi, for fixed 
B2 = 0.45. In all cases A = 1.0, C = 0.3, and E = Q. 
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0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 
temperature(T) temperature(T) 
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0.2 0.4 0.6 
temperature(T) 



0.4 0.6 O.E 
temperature(T) 



FIG. 6: Bend as a function of temperature for the system 
under an electric field, (a, b) For different values of the bend- 
polarization coupling C. (c, d) For different values of the 
electric field E. Parts (a) and (c) show systems with a biax- 
ial nematic phase, with parameters as in Fig. [3^, while (b) 
and (d) show systems with a direct transition from uniaxial 
nematic to polar, with parameters as in Fig. ^p. 



and P by averaging over the entire system. In the polar 
phase, because there is a systematic bend in the molecu- 
lar orientation, these global averages are reduced. That 
is the reason why Fig. [3] shows a decrease in S in the 
polar phase. Presumably this issue could be avoided by 
calculating local averages of the order parameters, but we 
use global averages because of the limited system size. 

To explore the phase diagram of the system, we carry 
out a series of simulations where we vary the coupling 
coefficients Bi and B2 and measure the resulting order 
parameters. Figure [5] shows two cross-sections through 
the phase diagram. Note that increasing B2 enhances 
the stability of the biaxial nematic phase at the expense 
of the uniaxial nematic and isotropic phases, but does 
not affect the polar phase. By comparison, increasing Bi 
enhances the stability of the polar phase at the expense 
of all the other phases. 
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FIG. 7: Equilibrium configuration from a Monte Carlo sim- 
ulation just above the uniaxial nematic-polar transition tem- 
perature. 



We now simulate the converse flexoelectric effect by 
applying an electric field E along the Z-axis. As ex- 
pected, the dipole direction b aligns parallel to the field, 
while the director n bends across the system, as a func- 
tion of X or F. We calculate this induced bend from 
Eq. (ITSI) as a function of temperature, electric field, and 
coupling coefficients. The results are shown in Fig. |6l 
Note that the bend increases as the field increases, as 
the bend-polarization coupling C increases, and as the 
temperature decreases. In particular, as the temperature 
decreases toward the transition into the polar phase, the 
bend responds sensitively to any applied field. The tem- 
perature dependence of the bend is sharpest for very low 
fields, and it is rounded off for larger fields. This be- 
havior is similar regardless of whether the system passes 
through the biaxial nematic phase or goes directly from 
uniaxial nematic to polar. These trends are the normal 
behavior of an induced order parameter in the presence 
of a symmetry-breaking field above a second-order phase 
transition. They are consistent with the Landau theory 
presented in the previous section, which suggests that 
the flexoelectric and converse flexoelectric effects become 
large near an incipient polar phase, where the uniform ne- 
matic phase is almost unstable and the system is most 
sensitive to any symmetry-breaking perturbation. 

To visualize the molecular configuration responsible for 
the large flexoelectric effect, we perform a simulation in 
the uniaxial nematic phase in zero field at T — 0.42, 
slightly above the transition into the polar phase at 
Tjjp K, 0.32. A snapshot of the configuration is shown 
in Fig. [T] Clearly there are local correlated regions with 
bend and polar order, but these regions do not extend 
across the whole system. At this temperature, the system 
is highly susceptible to an applied electric field, which or- 
ders the local correlated regions and induces a global po- 
larization and bend. Presumably it would also be highly 
susceptible to an applied torque coupling to the bend, 
which would have the same effects, although we have not 
done that simulation explicitly. 



Apart from the Monte Carlo simulations, the lattice 
model of Eq. p7)) can also be solved analytically through 
mean-field theory. The mean-field calculation is pre- 
sented in the Appendix, and the results are similar to 
the simulations presented in this section. 



IV. MODULATED PHASES: LATTICE MODEL 

As noted in the previous section, when the system 
cools into a polar phase with spontaneous polar order, 
it also acquires spontaneous bend. A simple example 
of the spontaneous bend is illustrated in Fig. [2Jl, which 
shows a gradual bend across the system, between the 
free boundaries. However, this simple configuration can- 
not be extended to give the molecular orientation across 
a larger system. In general, it is impossible to fill space 
with pure uniform bend. Rather, the system must form 
a more complex modulated phase, which might have a 
regular array of defect walls, or might have a mixture 
of bend with splay or twist in the director. Indeed the 
problem of filling space with bend in a polar liquid-crystal 
phase is quite analogous to the problem of filling space 
with twist in a chiral liquid-crystal phase, as discussed in 
Ref. [Hi. 

To determine the modulated structure of the polar 
phase, we repeat the simulations of the previous section 
with three modifications: we use periodic boundary con- 
ditions in all three directions, we use a slightly larger 
lattice of size 20 x 20 x 20, and we increase the bend- 
polarization coupling coefficient C to increase the bend, 
i. e. reduce the wavelength of the modulated structure, 
so that a full wavelength will fit in the simulation box. 
In these simulations no electric field is applied, so the 
only polarization is spontaneous order. We begin the 
simulations in the high-temperature isotropic phase and 
gradually cool into the low-temperature polar phase. In 
this way, the system is free to select its own modulated 
structure. 

In these simulations, two distinct types of modulated 
structures form, depending on the model parameters. 
The first structure, shown in Fig. [SJ is equivalent to the 
twist-bend phase proposed by Dozov Q. In this struc- 
ture, the director n has a helical modulation, which is 
randomly right- or left-handed. The director is not per- 
pendicular to the helical axis, as in a cholesteric liquid 
crystal. Rather, it precesses in a cone about the helical 
axis, with a fixed cone angle between 0° and 90°. As 
a result, the director deformation is a mixture of twist 
and bend, unlike a cholesteric liquid crystal which has 
pure twist. The dipole direction b also precesses about 
the helical axis, while remaining perpendicular to n and 
perpendicular to the helical axis. This structure is spa- 
tially homogeneous, in that every position is equivalent 
to every other position with a twist. Hence, every posi- 
tion has the same energy, and there are no defects. 

The second structure, shown in Fig. [SI is equivalent to 
the splay-bend phase proposed by Dozov 7]. Here, the 
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FIG. 8: Equilibrium configuration in a Monte Carlo simula- 
tion of the twist-bend phase, for model parameters A = 2, 
Bi = 0.5, B2 = 0.4, C = -2.0, and T = 0.5. (a) Top view, 
(b) Side view. 





FIG. 9: Equilibrium configuration in a Monte Carlo simula- 
tion of the splay-bend phase, for model parameters A = 3, 
Bi = 0.15, B2 = 0.4, C = -2.0, and T = 0.3. The director 
and polar order are both in the plane of the figure. The labels 
indicate the the local polarization along the vertical axis. 



and direction. In the pure-bend regions the b vectors are 
very well aligned, and hence the polar order parameter 
has a large magnitude. These regions are indicated by 
the labels P > and P < in the figure. By contrast, 
in the zero-bend regions, the b vectors are disordered, 
and hence the polar order parameter averages to zero. 
These regions are indicated by the label P = in the 
figure. The zero-bend regions have a higher energy than 
the pure-bend regions, so they can be regarded as defect 
walls. 

The simple example of Fig. [2ji can be understood as 
one pure-bend region going across the finite simulation 
cell. The splay-bend structure of Fig. [3] shows how this 
structure can fill space with a periodic alternation of 
pure-bend "defect-free" regions and zero-bend "defect 
walls." 

When we say that the structures of Figs. [S] and [S] are 
equivalent to Dozov's twist-bend and splay-bend phases, 
we mean that the director modulations are the same as 
what he proposed. His paper does not explicitly con- 
sider the polarization modulation, although his sketches 
suggest a variation in the polarization direction that is 
similar to our simulation results. The twist-bend phase 
has also been visualized in a simulation by Memmer f§\ . 



V. MODULATED PHASES: LANDAU THEORY 

To understand the modulated phases better, we return 
to the Landau theory of Sec. II. These phases occur for 
temperatures below the critical temperature of Eq. ((5]), 
where the quadratic form in the free energy of Eq. ^ is 
not positive-definite. Hence, we must add further terms 
to stabilize the free energy. First, there must be a term of 
ji/F'^, which keeps the magnitude of the polar order from 
increasing without limit. Second, there must be a term 
of ^k(VP)^, which penalizes gradients in the magnitude 
and direction of polar order. (In terms of tensor indices, 
we interpret this gradient term as ■^K{diPj){diPj). There 
could be other tensor contractions, but they do not mat- 
ter for our analysis.) With these new terms, the Landau 
free energy becomes 



F 



AB P 



+iA*'(T-r„)p2 + i^p4 + i«,(VP)2. (20) 

To model the twist-bend phase, we make the varia- 
tional ansatz for n{x) and P(a;) inspired by the simula- 
tion results of Fig. [H 



director h goes back and forth within the plane of the 
figure. As a result, the director deformation is a mixture 
of splay and bend. Note that different regions are not 
equivalent to each other — some regions have almost pure 
bend, and other regions have almost zero bend. In this 
structure, the local polar order varies in both magnitude 



= (l-a2)i/2 
P(.t) = — pcos((7a;)y +psin(ga;)z 



h{x) = (1 - a^)^''^x + asin(5a::)y 



! cos(ga;)2 



(21) 



This ansatz has three variational parameters: a is the 
sine of the cone angle for the director, p is the magni- 
tude of the local polar order, and q is the wavevector of 
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the modulation. In terms of these parameters, the splay, 
twist, and bend become 



S(x) = n(V • fi) = 0, 

T(x) = n- (V X n) = a^q, 

B(x) = n X (V X n) 

— a(l - a'^y^^q[— cos((7.T)y + sm{qx)2 



(22) 



Note that the splay is zero, as it should be for the twist- 
bend deformation. The twist and the bend magnitude are 
constants, while the bend direction precesses in a helix. 
Plugging those quantities into the free energy gives 

Ftb - ^K2aS^ + ^K3a\l-a^)q^ -Xapq{l-ay^^ 



, ^ I frr Ti\2|l 4,1 22 

+ [T - Tq)p + -vp + -up q . 



(23) 



Minimizing the free energy over the variational parame- 
ters a, p, and q then gives the behavior near the transi- 
tion, for T <Tc, 



P 



2A \K2 
K'ifi' ( 3 



4A2 \2K2n 



[Tc - T), 



(24) 



Also, the free energy of the twist-bend phase just below 
the transition is 



(25) 



By comparison, to model the splay-bend phase, we 
make the variational ansatz inspired by the simulation 
results of Fig. El 



n(a;) — C0S(/)(a;)x -|- sin(/)(x)z. 



(26) 



P(a;) — — pcos(ga;) sin 20(x)x — pcos((7a;) cos (/i(a;)z, 

where 0(x) = Q!sin(ga;). This ansatz also has three vari- 
ational parameters: a is the amplitude of the director 
modulation, p is the amplitude of the polarization mod- 
ulation, and q is the modulation wavevector. For this 
state, the splay, twist, and bend become 

S(a;) = — -(7acos(gx) sin2(^(a;)x — gacos(qa;) sin0(x)z. 



T{x) = 0, 



(27) 



B(a;) — — qa cos(qx) sin 2(/)(a;)x — go! cos((j'a;) cos 0(a;)z. 

Here the twist is zero, as it should be for the splay-bend 
deformation. The splay and bend both vary periodically 
through the modulated structure. Plugging these quan- 
tities into Eq. ([^0)1 gives a free energy density that also 



varies periodically through the modulated structure. We 
average the free energy density over the full modulation, 
and then minimize the average free energy over the vari- 
ational parameters a, p, and q. This minimization gives 
the behavior near the transition. 



A \Ki 
8A2 \2KiK 



(28) 



7 f^J^ 

2k 



1/2 



{T,-Tf'\ 



Furthermore, the free energy of the splay-bend phase just 
below the transition is 



FsB — — 



32XikA4 



(Tc - T)3 



(29) 



From these results, we see that the uniform nematic 
phase can become unstable to the formation of either 
the twist-bend phase or the splay-bend phase at the crit- 
ical temperature Tc- We can then ask which of these 
modulated phases is more stable. Comparison of the free 
energies ([^5)) and (UHl shows that the twist-bend phase 
is more stable if Ki > 2K2, while the splay-bend phase 
is more stable if Ki < 2K2 ■ Interestingly, this is exactly 
the same criterion for the relative stability of the phases 
calculated by Dozov Q • This criterion is reasonable, be- 
cause the elastic constants Ki and K2 give the energetic 
costs of splay and twist deformations, which are required 
in the splay-bend and twist-bend phases, respectively. 



VI. DISCUSSION 

In this paper, we have presented a theory for orien- 
tational order in bent-core liquid crystals. The theory 
combines three parts: the Oseen-Frank free energy for 
director gradients, a tendency toward polar order per- 
pendicular to the director, and a coupling between polar 
order and director bend. In Landau theory, the Oseen- 
Frank free energy is represented by the K terms, the 
tendency toward polar order by the ^ term, and the cou- 
pling by the A term. In the lattice model, the Oseen- 
Frank free energy is represented by A, the tendency to- 
ward polar order by Bi, and the coupling by C. Either 
way, the physical conclusions are the same: In the uni- 
form nematic phase, there is a flexoelectric effect, where 
an imposed bend leads to a polarization, and conversely 
an applied electric field leads to a bend. This flexoelec- 
tric effect increases as the temperature decreases toward 
a polar phase. At the critical point, the response to an 
applied field diverges, and there is a second-order tran- 
sition from the uniform nematic phase into a modulated 
polar phase. The modulation may have the twist-bend 
or splay-bend structure, depending on the relative values 
of Ki and K2- 
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To compare our theory with the previous work by Do- 
zov Q, note that two phenomena occur at the critical 
temperature Tc- the system transitions from the uniform 
nematic phase to the polar phase, and the renormalized 
bend elastic constant of Eq. (fTO|) changes sign from 
Kf > for T > Ta to Kf < for T < T^. Do- 
zov would say that the modulated phase is caused by the 
negative elastic constant K^^ . By contrast, we would say 
that the modulated phase and the negative K^^ are both 
caused by the bend-polarization coupling together with 
the tendency toward polar order. Of course, there is no 
contradiction between these two theories; they are just 
different ways of expressing the same physical concept. 

Our theoretical results can be compared with recent 
experiments. For the flexoelectric effect, the most rel- 
evant comparison is with the experiments of Harden et 
at. 0, , which found a surprisingly large bend flexoelec- 
tric coefhcient in bent-core liquid crystals, about three 
orders of magnitude larger than the typical value in rod- 
like liquid crystals. This observation is at least qualita- 
tively consistent with our concept that bent-core liquid 
crystals are near an incipient polar phase, and hence are 
very sensitive to any slight polar perturbations. How- 
ever, one aspect of this comparison is confusing. In our 
theory, it is easy to understand why the polarization and 
bend induced by an applied electric field or an applied 
torque should be very large in bent-core liquid crystals, 
because these quantities diverge at the critical point. By 
comparison, it is not easy to understand why the ratio 
between induced polarization and induced bend should 
be very large, because this ratio does not diverge at the 
critical point. (The ratio increases as T —> Tc, but only 
toward a finite limit.) Even so, that ratio is the standard 
definition of the flexoelectric coefficient, which Harden et 
al. found to be surprisingly large. One possible expla- 
nation of this discrepancy might be that the experiment 
involves local smectic order, which is not included in the 
theory; perhaps this smectic order increases the ratio of 
polarization to bend. An alternative explanation might 
be that the experiment is somehow measuring one of the 
response coefficients that diverges at the critical point 
and is not actually measuring the ratio of polarization to 
bend. Yet another possibility might be that the experi- 
ment is actually in a modulated polar phase that has not 
yet been identified. 

For the modulated polar phases, the most relevant 
comparison is with recent experiments that find nonuni- 
form nematic phases in systems of bimesogens 0, [lo| . 
These experiments provide good evidence that the ob- 
served modulation is a twist-bend structure, which is an 
encouraging consistency between theory and experiment. 
However, one point is confusing in the comparison with 
Ref. [^]. The experiment reports that the periodicity of 
the observed stripe pattern is twice the thickness of the 
cell. By contrast, the theory predicts that the periodic- 
ity should be determined by material parameters, even 
in the bulk liquid crystal, regardless of the cell thick- 
ness. One possible explanation of this discrepancy might 



be that surface anchoring modifies the predictions of the 
theory and fixes the periodicity of the incipient modu- 
lation. Such surface effects could be a topic for future 
research. 

As a final point, we note that the modulated polar 
phases are locally ferroelectric; they have spontaneous 
polar order which leads to electrostatic polarization. This 
local polarization is modulated in a helix for the twist- 
bend phase or a planar wave for the splay-bend phase, 
and hence it averages to zero globally. For that reason, 
it might be difficult to observe in an experiment. In this 
respect, the modulated polar phases are similar to ferro- 
electric smectic-C* liquid crystals, which also have local 
polar order that is modulated in a helix and averages to 
zero globally. By analogy with ferroelectric smectic-C* 
liquid crystals, there might be ways to unwind the helix 
of the twist-bend phase (or eliminate the wave modula- 
tion in a splay-bend phase) to obtain a structure with 
long-range polar order. For example, strong surface an- 
choring might give a surface-stabilized ferroelectric ne- 
matic phase. This surface stabilization could be another 
topic for future research. 
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Appendix: Mean-field calculation for lattice model 

The purpose of this Appendix is to show how the lat- 
tice model of Sec. HI can be solved analytically through 
mean-field theory, rather than Monte Carlo simulations. 
The calculation is analogous to the mean-field theory in 
our previous paper [J , but now for bend instead of splay 
flexoelectricity. The results are similar to the simulations 
presented in Sec. HI. 

For the mean-field calculation, we consider a small re- 
gion of the system in the uniform nematic phase under an 
applied electric field. We suppose that the field and av- 
erage polarization are along the x axis, and the director 
has a bend from layer to layer in the lattice as function 
of z. Hence, we write 

{(— sin A^, 0, cos AS) in layer 1, 
(0,0,1) in layer 2, (A.l) 

(sin A6', 0, cos A6') in layer 3. 
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FIG. 10: Mean-field results for the order parameters showing 
different types of transitions, (a, b) The biaxial order param- 
eter P2 (•) and polar order parameter Pi (o) as functions of 
temperature T, for Bi = 0.03 (in a) and Bi = 0.20 (in b) at 
zero field, (c, d) Bend Ad (*) as a function of temperature 
under an applied field. In all cases A — 1.0, B2 = 0.45, and 
C = 0.30. 



Similarly, for the dipole direction on site i, we write 



(cos AO cos 6i , sin 9i , sin cos 6i ) 
b,: = ( {cos9i,sm9i,0) 

(cos A9 cos 9i, sin 9i, — sin A9 cos 9i) 



in layer 1, 
in layer 2, 
in layer 3. 
.(A.2) 

We now suppose that the uniaxial nematic order is per- 
fect, with a single well-defined value of the bend angle 
A0, but there is a statistical distribution of the dipole 
orientations 9. In terms of this distribution, the polar 
order parameters is Pi = {cos9) and the biaxial order 
parameter is P2 = {cos29). If we average the lattice 
Hamiltonian of Eq. p7)) over the statistical distribution 
of 9, and assume that A9 is small, then we obtain the 
average energy per site 



U = -A[3- {A9)^] - BiPl 



■i-\{A9f 



-Bo 



-{1 + P2)\A9)' 



(A.3) 



-CPiA9 - EPi 



To construct the free energy, we must combine the av- 
erage energy with the entropy associated with the distri- 
bution of dipole orientations 9. The orientational distri- 
bution function is 



p{9) = 



^(v\ COS 9-\-V2 COS 29) 



J^'^ g(t;i cos e-\-V2 cos 2e)fj^0 



(A.4) 



where vi and V2 are two parameters in the effective poten- 
tial acting on 9. They are related to the order parameters 
Pi and P2 by 



Pn{vi,V2) = 

for n = 1 and 2. 
then becomes 



/g cos(n^)e("i ^+''^ """^ 



J^^ g(«i cos6»-l-U2 COS 261) 



(A.5) 



The entropic part of the free energy 



2ir 



TS = keT / pi9)\og{p{9))d9 



ViPi + V2P2 



(A.6) 



log 



2tt 



^{vi cos 9+V2 COS 26)^^ 



Combining Eqs. (jA.3[) and (jA.6p . the full free energy is 
F = U -TS. 

The problem to be solved is now: For a given set of 
interaction parameters A, Bi, B2, C, E, and temperature 
T, we should find the variational parameters Vi, V2, and 
A9 that minimize the free energy. Once we know the 
values of vi and V2 , we can calculate the order parameters 
Pi and P2. We can then determine whether the phase 
is uniaxial nematic (Nu), biaxial nematic (Nb), or polar 
(P), and we can see how the bend A9 is related to the 
order parameters. 

The free energy is minimized with Mathematica. Fig- 
ure [TOb.b shows plots of the order parameters and the 
bend in zero field as functions of temperature. Clearly 
the qualitative behavior of the phase transitions found 
in Monte Carlo simulations is reproduced. For small Pi 
and P2, there is a transition from Nu Nb, followed 
by a transition from Nb — P at lower temperature. For 
slightly larger Pi , the two transitions merge into a direct 
transition from Njj — > P. These mean-field results are 
generally consistent with the simulation results shown in 
Fig. m except that the transition temperatures are sig- 
nificantly higher. As mean-field theory exaggerates the 
tendency towards an ordered phase, it overestimates the 
transition temperatures. Despite this limitation, the sim- 
ple mean-field calculation is successful in describing the 
qualitative behavior of the system. 

To see the effect of an applied electric field, we repeat 
the minimization with a small field E in the free energy. 
FigurefTOh.d shows the variation of the bend as a function 
of temperature under a field. These mean- field results are 
generally consistent with the simulation results shown in 
Fig. [6l showing a bend that responds more sensitively 
to any applied field as the temperature decreases toward 
the transition into the polar phase. 
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